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Abstract 

We consider the spin orientation of the final Z bosons for the processes in the 
Standard Model. We demonstrate that at the threshold energies of these processes 
the analytical expressions for the Z boson polarization vectors and alignment tensors 
coincide (e + e~ — > ZH, Zj) or are very similar (e + e~ —> ZZ). In addition, we present 
interesting symmetry properties for the spin orientation parameters. 



1 Introduction 



The Standard Model (SM) of elementary particle physics has been phenomenologically 
very successful. However, as a theory it is less impressive. There are many fundamental 
questions that remain unanswered by the SM. Due to this, there are strong reasons to 
expect that there is a great deal of (new) physics beyond the SM, more fundamental with 
a characteristic mass scale. The present good agreement between the accelerator-based 
experimental data and the SM predictions suggests that the energy scale associated with 
any new physics model should be much higher than the energy reach of the colliders up 
to now. It is believed that the new accelerator, the Large Hadron Collider (LHC), started 
in 2009, is able to produce new physics particles, first of all new on-shell resonances or a 
single heavy new particle in association with the SM ones [TJ. However, such a scenario 
for visualizing new physics experimentally is obviously beyond the reach of future e + e~ 
colliders, like the planned International Linear Collider (ILC). 

There are methods to probe new physics at energies below the new physics mass scale. 
These more indirect scenarios are based on observations of small deviations from the SM 
predictions for processes with SM particles where new physics effects can arise only from 
non-standard interactions. The price to pay for such a possibility is the need for higher 
sensitivity, both experimentally and theoretically. In the case of future ILC-type colliders, 
the achievable sensitivities are expected to be dramatically (i.e. up to orders of magni- 
tude) improved as compared to the ones obtained up to now. Among the various ways 
to increase the sensitivity we emphasize two. First, there is a reason to expect that pro- 
cesses including heavy particles are more sensitive to new physics manifestations. Second, 
processes including polarization effects have higher sensitivity than unpolarized processes. 
The possibility at ILC to use both longitudinally and transversely polarized initial beams 
provide additional means to increase the sensitivity for new physics. 

Aforementioned circumstances considerably increases the potential to probe anomalous 



2 



interactions via spin effects in e + e~ annihilation into heavy particles. Therefore, the pro- 
cesses e + e~ — > ZH, Zj, and ZZ have been studied extensively. At the same time it is 
obvious that the search for new physics, even though performed at very sensitive future 
colliders, cannot be successful without knowing the SM predictions with appropriate pre- 
cision. In Ref. pQ Lykken writes: "[. . .] to first approximation [. ..] experimenters do not 
need to know anything about BSM [Beyond the SM] models in order to make discoveries 
- but they need to know a lot of the Standard Model physics!" In the context of search- 
ing for new physics in aforementioned processes via spin orientation, this means that the 
knowledge of all possible spin effects in these processes in the framework of the SM forms 
a good basis for rejecting or limiting anomalous couplings. In this paper we try to add our 
contribution to this basis. 

In the framework of the SM we present and analyse analytical expressions for the Z 
boson spin polarization vectors and alignment tensors for the processes e + e~ — > ZH, Zj, 
and ZZ near the production threshold. For all these processes the initial beams are taken 
to have both longitudinal and transverse polarization components. We demonstrate that 
the expressions for the spin orientation parameters in these processes have similar shape 
and symmetry properties. Even more, the expressions for the Z boson polarization vector 
and alignment tensor for the processes e + e~ — > Z'j and ZH coincide. 

The paper is organized as follows. In Sec. 2 we explain how spin orientation parameters 
can be described by density matrices. In addition we show how to obtain the actual spin 
polarization vector and alignment tensor of the final Z boson. In Sec. 3 we present the 
analytical expressions for the spin orientation parameters of the three processes and dwell 
on similarities and interesting symmetry features. The results are checked with the help 
of positivity conditions. At the end of this paper we give our conclusions. 
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2 Description of spin orientation parameters 

In physical processes the particles with non-zero spins are generally in a mixed spin state. 
Contrary to the pure state which can be described by a single wave function, the description 
of a mixed state deserves an incoherent superposition of 2s + 1 orthogonal pure state wave 
functions. By definition this superposition means that in order to calculate the probability 
for finding a certain mixed state, one has to calculate the average of the probabilities of 
the pure states, assigning to each pure state an appropriate weight. 

The most natural way to describe mixed states is to use the spin density matrix formal- 
ism. Technically, if interested in the spin orientation of a certain particle in the process, 
one replaces the element uu (spin-1/2 case) or e^e* u (spin-1 case) in the squared amplitude 
of the process by the mixed state spin density matrix. For our calculations the relativis- 
tic spin-1/2 density matrices for the initial beams and the non-relativistic spin-1 density 
matrices for the final Z-bosons are needed. The relativistic spin-1/2 density matrices have 
the familiar shape 



where the upper signs refer to the electron and the lower ones to the positron. Here and 
in the following we use Feynman's slash notation $ = . The polarization four- vectors 
are given by 



where £ T are the polarization vectors in the resp. rest frame of the particle. For the absolute 
value |£ T | of the mixed state polarization vectors one can take any value between zero and 
one. If the electron and positron have both longitudinal polarization (LP) and transverse 
polarization (TP), it is useful to divide the polarization vector s T up into LP and TP parts. 
The limit m/ko — > which is used in our calculations can be conveniently taken by making 
use of the approximation 



Pt = 2^=f ±"0(1 +75^ T ) 



(1) 






m 



(3) 
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and subsequently setting m = 0. h T is the measure of LP and r T = (0, r T ) is the TP 
four- vector with f T the transverse component of the polarization vector (k T ■ r T = 0). 
After inserting Eq. (J3j) into Eq. (CTJ and performing the limit m — > 0, one obtains [2] 

P = ^(l±M5 + 7 5 fJ T (4) 

convenient for our calculations. 

For describing mixed states by 3 x 3 non-relativistic hermitian spin-density matrixes, 
at most 8 real parameters are needed. For expanding such a matrix one can use a basis 
containing the unit matrix /, the three spin matrices Si, and the six combinations of 
products of two spin matrices in the form 



Sij — — ySiSj + SjSi — fiy ; 1 . (5) 

In such a basis the basic elements are hermitian zero-trace matrices which makes their 
usage convenient. The density matrix expanded in this basis reads 



p=^(l + h i S i + h ij S ij ) (6) 



where here and in the following we sum over all indices which occur twice (Einstein's 
convention). Note that this basis is over-determined. Instead of 9 matrices needed to 
expand the density matrix it contains 10 matrices. However, not all the elements of the 
basis are linearly independent. The relation 

S X X ~\~ Syy + S ZZ = (7) 

holds for three of the basic elements. In order to reduce the number of parameters from 9 
to 8 independent ones, it is necessary to claim 

txx tyy ~\~ t zz 0. (8) 



Due to this constraint, one obtains 



Pij — -^ij + -tk(Sk)ij + ~tki(SkSi + SiSk)ij (9) 



where tk = Tr(pSk) is the polarization vector and %i = Tr(pSki) is the orientation or 
alignment tensor describing the alignment of spins. 

When taking the spin-1 matrices in the representation (Sk)ij = —itijk where is the 
totally antisymmetric tensor, one obtains 

(SkSi + SiSkjij = 25ki5ij — 5ki5ij — 5tj5u. (10) 

In this case, pij from Eq. ([9]) can be written in the simple form 

Pij = 3 \ - -jitk^ijk -Ujj ■ (11) 

This form is used in our calculations. 



2.1 Spin orientation of final particles 

In our calculations it is supposed that the spin orientation of a single final Z boson is 
observed!] Accordingly, the squared amplitude of the process can be given in the form 

\M\ 2 -S + ViU + TijUj (12) 

where S, Vi and are the scalar, vector and tensor built from the polarization param- 
eters of the initial electron and positron, the kinematical parameters of all the particles 
participating in the process, and the coupling constants. Since the calculations are per- 
formed in the CM frame with vanishing electron mass at threshold energies, there are only 
two kinematical parameters - the electron three-momentum k and the mass Mz of the Z 
boson. 

lr The spin orientation parameters U and Uj determine the angular distribution of the decay products. 
Therefore, the Z boson spin orientation parameters can be seen by analyzing the angular distribution of 
its decay products. According to the Particle Data Group [8], the Z boson decays with 10% probability 
into lepton pairs (e + e _ , fj, + fi~, t + t _ ), and with 70% probability into hadron pairs. Considering only 
lepton pairs which can be clearly distinguished from the background, the analysis can be found in Sec. 4 
of Ref. [3] , giving the correlation between the spin orientation parameters and the angular distribution. 
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The squared amplitude determines the probability that the process produces a final Z 
boson with spin orientation parameters ti and ty. On the other hand, the same probability 
can be expressed also as the traced product of two density matrices 

Tr(p» ~ 1 + hlU + (13) 

where the real (actual) density matrices and their parameters are denoted by an upper 
index r. The symbols without this index are the density matrices and their parameters 
that are substituted into the squared amplitude instead of s^e^* . By comparing Eqs. ( 112"]) 
and ffl3|) one can find the actual polarization vector and alignment tensor of the Z boson, 

** = ^> $ = fCty + ra<,-) (14) 

where the term T incorporates the tracelessness (JSj). 
2.2 Positivity conditions 

The values of the parameters ti and tij are restricted by their physical boundaries as 
extremal values. In this physical region parameters are linear independent. However, due 
to the positivity (non-negativity) condition of the density matrices they depend on each 
other non-linearly. It can be shown that the most restrictive requirements for the spin-1 
density matrix can be given by the inequality [4] 

2Trp 3 - 3Trp 2 + 1 > 0. (15) 

If one expresses the density matrix through its parameters, one obtains 

2 1112 

9 _ 2^- 7 ^- 7 _ g^i^J — ^tijtjktki > 0. (16) 

From this expression one can clearly see how the vector and tensor parameters depend 
non-linearly on each other. Eq. ( fl6|) is useful when analyzing the possibilities of tuning the 
final Z boson spin orientation by varying the polarization of the initial beams. 



3 Z boson spin orientation parameters at threshold 

We are now ready to present the final Z boson spin orientation parameters (polarization 
vectors and alignment tensors) at the threshold energies for the three processes e + e~ — > 
ZH, Zj, ZZ, calculated in the center-of-mass system. Note that the kinematics simplifies 
significantly when calculating at the threshold. On the other hand, the results obtained 
are also valid close to the threshold. 

The first process e + e~ — > ZH could be important in clarifying the mechanisms of 
electroweak symmetry breaking proposed by various beyond-the-SM models (for an inde- 
pendent work on this subject see e.g. Refs. E]). In the SM at tree level this process 
is described by an s-channel Feynman diagram with a point-like ZHZ vertex (Fig. [1]). 
Actually, at lowest order there are three Feynman diagrams corresponding to the process. 




Figure 1: s-channel SM tree level diagram for e + e — > ZH with point-like ZHZ vertex 

However, the two diagrams where the Higgs boson couples to the electron line can be 
ignored due to the smallness of the coupling which is proportional to the electron mass. 

The two remaining processes e + e~ — > Zj, ZZ have been most of all used to investigate 
possible new physics manifestations via anomalous gauge boson self- couplings Z77, Z^Z, 
and ZZZ. In the SM at tree level these processes are described by a t-channel and an 
it-channel diagram (Fig. [2]). 

The Z boson spin polarization vectors and alignment tensors for the processes e + e~ — > 
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k + — * P2 — - k + 
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Figure 2: (a) t-channel and (b) w-channel SM tree level diagram for e + e — > Z7, 

Zif , Z7 with vanishing electron mass have been calculated by us earlier [31 II]- At the 

threshold energies the spin orientation parameters reduce to the same expressions, 
; ( 2 ). 



Hhres 



(2) 



k. 



t 



ij thres 



hi.i. _ is- ) + 39l9r 

2^ K 3 % 3) + ^.(2) 



1 t_ ■ r_ 



(17) 
(18) 



where we used 9l = ^(9v + 9a), 9r = \{.9v~9a) in order to write the expressions in a more 
symmetric form, and the short hand notation f± := <?f,(l— /t_)(l+/t+)±<7^(l+/i_)(l— 
/i T and r T are the measures of the electron (positron) longitudinal polarizations and the 
electron (positron) transverse polarization vectors, respectively, k indicates the unit vector 
along the electron momentum. 

Even though somewhat peculiar, the coincidence of the Z boson spin orientation pa- 
rameters of a s-channel diagram with the spin orientation parameters of a t-channel and 
an w-channel diagram at corresponding thresholds is expected. One can say that Eqs. ( IT71) 
and ( lisp give the spin orientation parameters of the real Z boson in the process e + e~ — > Z. 
The Higgs boson is a spin-0 particle and, therefore, cannot affect the spin orientation of 
the Z boson at the threshold of e + e~ — > ZH. In the same way, for pi — > the photon in 
the process e + e~ — > Z7 can be considered as a radiative correction to the main process 
e + e~ — > Z. The coincidence of the Z boson orientation parameters for the two different 
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processes, therefore, can be taken as a cross-check for the final expressions. 

In order to obtain the spin orientation parameters for one of the final Z bosons in 
the process e + e~ — > ZZ at threshold, we again take the (massless) electron and positron 
beams to be polarized, having both longitudinal and transverse polarization components. 
In summing over the spin orientations of the second Z boson (in Fig. [2] indicated with 
momentum P2), we obtain for the squared amplitude of the process at threshold 

3sm [2Vw) I 4 4 

+ 9l9r 



3sin 2 (2^) I ^ ^ 

(f_ ■ f + )kikj + T_iT +j + T-jT+i tjjj (19) 



where the short-hand notation f± := g\{\ — h_)(l + h + ) ± g%(l + h_)(l — h+) is used. 



According to the rules given in Sec. 12.11 for the spin orientation parameters one obtains 

/ (4) - 

^thrcs = Z [.(4) ^' (^^) 



•■> 2 J. 



(21^ 



Note that the terms 5^- are introduced into in order to make it traceless. Obviously, 
the spin orientation parameters in Eqs. (120]) and fT2"T]) do not coincide with those of the 
processes e + e~ — > ZH^Z^ (Eqs. (|17[) and (|18|) ) but are similar to them. This is quite 
understandable. Due to the two Z boson final state the constants gi and gn appear twice 
in the squared amplitude. Therefore, the squares in /± had to be replaced by fourth 
powers in f±K On the other hand, only the direction k of the electron and the transverse 
polarization vectors t± are available to build the spin orientation parameters. As a result, 
besides the change of the power of the coupling constants the expressions can differ only 
by the constants in front of the various three contributions. This is indeed the case in the 
expressions obtained. 
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3.1 Symmetry properties of the polarization vectors 

In this section we demonstrate that the polarization vectors of the Z boson in the different 
processes under consideration feature interesting symmetries. Expressing the polarization 
vectors as a function of the effective polarization parameter 

< 22 > 

of the initial beams, for the processes e + e _ — > ZH, Z7 we obtain t t f res = t z (x)k where 



t z ( X ) = a = = 4^ « 0.147 (23) 

ax + 1 9l + 9r 9y + 9 a 



where we used gy = |(— 1 + 2 sin 2 6w) ~ —0.037, = — | = —0.5, gi = — § cos 2 9 W , and 
gn = \ sin 2 ^vf (with sin 2 = 0.2315, cos 2 0w = 0.7685 [8]). For the process e + e~ — > ZZ 
one obtains t t f^s = t ZZ (x)k where 

t zz (x) = -f^k 

b = = ( t 9V9 i 9 j + A 9 \\ - 0.294. (24) 

9t + 9% (g v + g A ) 2 + 4g v g A 

For longitudinally unpolarized initial beams (x = 0) one obtains quite similar values t z (x = 
0) = a ~ 0.147 and t zz (x = 0) = |6 ~ 0.144. The similarity is related to the smallness 
of gy. Exact coincidence would hold for gy = 0. However, in this case one would obtain 
a = b = 0. 

On the other hand, the polarizations vanish for % = —a and x = ~b, respectively, i.e. 
t z (x = ~ a ) = an d t zz (x = —b) = 0. This is because the reciprocal functions read 

S® = ±±i, *"« = g££ (25) 

Note that for the process e + e~ — > ti at threshold the SM result for the final top quark 
polarization can be expressed as £ = £( X )k where the function = — {x + C )/( C X + 1) 
has similar symmetry properties with £(x = 0) = x(£ = 0) = — c ~ —0.408 [9J. 
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3.2 Positivity condition for longitudinally polarized initial beams 

In Fig. [3] the dependence of t z and t zz as functions of the effective polarization parameter 
X is presented. Note that \ can obtain values between —1 and +1 only if the initial beams 
have no transverse polarization components. While t z takes all values between —1 and 
+1, the range for t zz is restricted to \t zz \ < |. This is mirrored also by the positivity 
condition in Sec. 12. 2\ giving a possibility to check our result in an independent way. 

As mentioned before, the polarization vectors and the alignment tensors can depend 
only on the direction k of the electron beam and the transverse polarization vectors of the 
initial beams. Considering first longitudinally polarized beams (t± = 0) and choosing the 
z-axis to be the direction of the electron beam, one obtains the non- vanishing components 

,z X + a i _ X + a 
ax + 1 " ax + 1 

Inserting the Uj into the positivity condition ( fl6|) it turns out that the condition is satisfied 
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independently of the value of t z and, therefore, also independently of the value of x- Hence, 
by varying the effective polarization parameter x one can theoretically force the Z boson 
polarization vector to take any value between —1 and +1. Practically achievable values are 
not far from the theoretical ones. When using the values h- = ±0.8, h + = ±0.6, planned 
to be achieved at the ILC, one can reach value x = ±0.95 for the effective polarization 
parameter, leading to similar values for the Z polarization (cf. Fig. EJ^a)). 

When doing the same with the components of the spin orientation parameters for the 
process e + e~ —> ZZ, one obtains 



az 1 X + b r _ 1 x + b 



2b X + l 2b X + V 

1 rZZ _ ,ZZ . £ 

2' xx yv ~ 4' 



J.ZZ j_Z Z j.Z Z {0*7} 

L zz ~ o' xx ~ L yy ~ ~- \ ZI ) 



Inserting the tij into the positivity condition ffTB 7 ]) . one obtains \t^ z \ < 1/2 which is consis- 
tent with Fig. [3]^b). 



3.3 Positivity condition as a test for the result 

Also in the general case, i.e. with transversely and/or longitudinally polarized initial beams, 
the positivity condition fflBT) is in harmony with the spin orientation parameters £j and 
of the processes e + e~ — > ZH, Zj, ZZ. Our investigation checks the obtained expressions 
for the spin orientation parameters and shows how much one can tune these parameters 
when varying the initial beam polarization parameters h T and f^. 

Taking both transverse polarization vectors r_ and f + to be directed along the x-axis 
of the coordinate system, for the processes e + e~ — > ZH, Z'y one obtains 



t z - 


f m 




/f 




t z - 

XX 


1 


3gL9R 


"2 + 


f (2) 
J + 


t z - 

yy 


1 


3gLgn 


_ 2 ~ 


/f 
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4 = 1- (28) 

Inserting the into the positivity condition f fT6|) . one finds that the condition takes the 
form = 0. This means that the positivity condition is satisfied in any case, i.e. for any 
values h T and t tx . As a consequence, by varying the LP parameters h T (or x) one can 
tune the Z boson polarization over the whole range [—1,1] if t_ or r + or both of them 
are equal to zero. In the same way, by varying r T one can tune the alignment tensor for 
h_ = h + = 0. 

Under the same conditions, for the process e + e~ — > ZZ one obtains 



t zz 



2 /f 



t zz _ 1 , M& T T 
L xx - 4 -I- 2 ^, (4 ) r -xr+x, 

.ZZ _ 1 ggjgL 

r ra ~ 4 2 j|4) r -^ r +^' 

4 Z = ~\ (29) 
After inserting £^ into the positivity condition (fTE]) . one obtains the restriction 

&) 2 < i - #(r-,r + ,) 2 . (30) 
/+ 

In absence of the term depending on t^ x t +x this result reproduces the previous result, 
i.e. that by varying the parameters h T (or x) one can theoretically tune the Z boson 
polarization from zero to 1/2 (t^ z = for x — ~ b = —0.288 and t zz = 1/2 for x — !)■ 
Finally, if we insert also the value t z zz into the positivity condition f fT6|) . we obtain a 
condition of the form 

(r- x r +x ) 2 < (l-/^)(l-/4). (31) 

For /i_ = h + = the condition (r_ 1 .r +x ) 2 < 1 puts no constraints on the transverse polar- 
ization of the initial beams. On the other hand, if one of the initial beams is longitudinally 
polarized, restrictions on the transverse polarizations of the initial beams are imposed. 
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4 Conclusions 

We have presented the spin orientation parameters of the final Z boson for the three 
processes e + e~ — > ZH, Z7 and ZZ in the Standard Model at the threshold energies 
and have found their properties. Note that a part of the properties found are valid also 
at higher energies. Therefore, the processes depend on the longitudinal polarization of 
the initial beams via the factors f±^ = g 2 L {l — h_)(l + h + ) ± g 2 R (l + hS)(l — h + ) or 
/± = g£(l — /i_)(l + /i + )±g^(l + /?,_)(l — /i + ) not only at the thresholds of the processes but 
also in the general case. Such properties could be helpful in separating SM contributions 
from new physics ones. 
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